Oscillation effects on neutrino decoupling in the early universe 



O 

o 

(N 



Steen Hannestad 
NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
(February 1, 2008) 

In the early universe, neutrinos decouple from equilibrium with the electromagnetic plasma at 
a temperature which is only slightly higher than the temperature where electrons and positrons 
annihilate. Therefore neutrinos to some extent share in the entropy transfer from e"'"e~ to other 
species, and their final temperature is slightly higher than the canonical value = [A/llY^^T-y. 
We study neutrino decoupling in the early universe with effects of neutrino oscillations included, 
and find that the change in neutrino energy density from e^e~ annihilations can be about 2-3% 
higher if oscillation are included. The primordial helium abundance can be changed by as much as 
1.5 X 10~* by neutrino oscillations. 



PACS numbers: 95.30.Cq,14.60.Pq,14.60.Lm,13.15.+g 
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Freeze-out of particle species from thermal equilibrium is one of the most important features of the early universe. 
One example is supersymmetric cold dark matter such as the neutralino. If these particles were kept in thermal 
equilibrium throughout the history of the universe their present day abundance would be suppressed by a huge factor 
because of their large mass. However, because of their weak interactions, the pair annihilation processes which keep 
the particles in equilibrium stop being efficient as a temperature of the order T ~ m/20 0]. At this time the particles 
decouple from equilibrium and their abundance is afterwards only diluted by cosmological expansion. Therefore such 
particles can be sufficiently abundant today to make up the dark matter 

In the present paper we discuss the decoupling of neutrinos from thermal equilibrium at a temperature of 1 
MeV. In the canonical picture the electron neutrinos decouple at a temperature of about 2 MeV, whereas the /xand r 
neutrinos, which do not have charged current interactions, decouple at roughly 4 MeV. Shortly after this decoupling 
the temperature drops below the rest mass of the electron and electrons and positrons pair annihilate and dump their 
entropy into photons. The neutrinos are not supposed to share in this entropy transfer and therefore their temperature 
is lower than the photon temperature by a factor (4/11)^/'^ [||. It is also implicit in this treatment that oscillations 
between active neutrino species have no effect because their distributions remain identical. 

However, because the temperature difference between neutrino decoupling and electron-positron annihilation (at 
T ^ Trie/ 3) is only one order of magnitude the neutrinos will to some extent share in the entropy transfer and the 
final temperature of the neutrinos will be slightly different from the canonical value. 

This problem has been treated many times in the past with various methods P-[T5[. Early calculations made 
various approximations , whereas recent treatments have included solving the full momentum dependent system 
of Boltzmann equations including all neutrino reactions [^13 



In all treatments so far, however, neutrinos have been assumed to be non-mixed species (although the possibility that 
neutrino heating could be changed by oscillations was mentioned in Ref. [^6)). The recent results from the Sudbury 
Neutrino Observatory (SNO) [0 experiment have confirmed that the Solar neutrino deficit is indeed explained by 
active-active neutrino oscillations, most likely — v^i oscillations JlSt . However, the specific values of the mixing 
parameters have not been conclusively measured. At present there are four different solution regions in mixing 
parameter space. The best fit points for these four solutions (from Ref. are listed in Table I. 

TABLE I. Best fit values of mixing parameters for solar neutrino solutions, as well as their goodness of fit. 



Solution 


(SmVeV^ 


sin 26o 


goodness of fit 


Large Mixing Angle (LMA) 


4.5 X 10"'' 


0.91 


59% 


Small Mixing Angle (SMA) 


4.5 X lO"'^ 


3.94 X 10"^ 


19% 


Low 


1.0 X lO"'' 


0.99 


45% 


Vacuum (VAC) 


4.6 X 10""' 


0.91 


42% 
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Also, by far the best explanation for the atmospheric neutrino problem is oscillations with maximal mixing between 
and another neutrino, most likely Vr, as suggested by the Super-Kamiokande experiment |Q. In the present 
paper we treat neutrino decoupling in the early universe, taking into account neutrino oscillations. Unfortunately 
the numerical complexity of the problem increases dramatically if oscillations are introduced. We therefore treat the 
problem in the so-called quantum rate equation approximation where the momentum dependence of the problem is 
integrated out. We also refrain from treating the full three-flavour oscillation problem and concentrate on Ve — 
oscillations while assuming that Vt is unmixed. Data on atmospheric neutrinos suggest that and Vt are in fact 
maximally mixed, and our assumption is therefore not likely to hold. Wc will discuss this issue more thoroughly in 
Section IV. 

We find that neutrino heating is slightly more efficient if oscillations are included. For non-oscillating neutrinos 

we find an increase in neutrino energy density of I^N^, ~ 0.0467, where AiV,y = 5p/ and p* — \ (rr)*^'^ P'y 
case of very efficient mixing of v^. and this is changed to l^N^, ~ 0.0479, a change of 2.5%. We also find that 
oscillations can change the primordial helium abundance. Neutrino heating induces a change of AYp ~ 1.2 x 10"'' for 
non-oscillating neutrinos, for maximal mixing this is changed to AYp ~ 2.2 x 10~^, i.e. a difference comparable to 
the total magnitude of the effect. In the limiting case of complete mixing of all neutrinos the result is ATV^ ~ 0.0486 
and AYp ~ 2.5 x 10"^. 

The next section contains a discussion of the essential rate equations needed to treat neutrino decoupling, as well as 
the equations needed to treat the cosmological expansion. Section III concerns the numerical details and results from 
solving the equations, and Section IV contains a discussion and conclusion. Finally, the Appendix contains various 
mathematical detail. 



II. BOLTZMANN EQUATIONS 



The fundamental equation describing the evolution of all particle species in the early universe is the Boltzmann 
equation, which for a non-mixed species is ||2l| ] 

i 

where CJ/] are the collision terms. In a homogeneous and isotropic system df /dt = df /dt. For the moment we do 
not consider the expansion of the universe, but this will be discussed at the end of the section. 

However, for a treatment of oscillating neutrinos it is necessary to replace the single particle distribution function 
with the density matrix, p, for the mixed species. For neutrinos the density matrix is a 3 x 3 matrix, but in the present 
work we limit ourselves to studying two-ffavor oscillations, with the third neutrino being non-mixed (in practice we 
study Ve — Vfj, oscillations with being non-mixed). The density matrix for neutrinos can then be written as 

p(p) = l[Po(p) + Pb)-a], (2) 

where ai are the Pauli matrices. An equivalent expression holds for anti-neutrinos. Notice that this definition is 
slightly different from the notation which is most often used p{p) = ^Pq{p)[1 + P{p) ■ a]. This notation has the 
advantage that P describes the internal state of the mixed neutrinos, whereas in our notation it does not. On the 
other hand our notation significantly simplifies the equations because it is linear (i.e. it does not contain a PqP term). 
The convention we use is the same as that used by McKellar and Thomson in their treatment of active-active 
oscillations in the early universe. The material in this section in follows their treatment closely in several respects. 

The usual one-particle distribution functions are the diagonal elements of the density matrix so that, na{p) — 
UMp) + Pzip)] and npip) = i[Po(p) - PAp)]- 

The mixed neutrinos, described by [Pq.P) then evolve according to the equations 

P{p) = V X P + [R^{p) - Rp{k)]i - D{p)Pt{p) (3) 
+ j dp'd{p,p')VT{p') - C(p)Po(p) + j dp'c{p,p')Po{p') (4) 
Po{p) = Rc.ip) + R0{p), (5) 
where Pt = Pe(p)x -|- Py(p)y is the transversal component of the polarization vector P. The rate terms are given by 
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Ri{p) = J dp'dkdk' 
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F,,{pk\p'k')[n,{k')nj{p') - n,{k)n,{p)] (6) 



-\Y^Gi{k'p'\kp)VT{p)-P*T{k) 



(7) 



where is over ah weakly interacting species and is over all other particles than the mixed neutrinos. In all 

cases / dq is taken to mean integration over phase space in the sense J dq — J ffj^i ■ 

Here, the first term in the brackets is the usual Boltzmann equation with exactly the same structure as for a non- 
mixed species. The second term arises due to the possibility that mixed state neutrinos annihilate with mixed state 
anti-neutrinos. This term therefore is only present for active-active oscillations where both states are interacting. The 
matrix element terms F and G are given by 

F,,{pk\p'k') = 2T:NV\]{p)J{k)\i{p'):i{k'))6{p + k-p' -k') (8) 
Gi{pk\p'k') = 2nNV{v^{p'),v^{k')\l{k)~l{p)) (9) 
xV{vp{p'), Dp{k')\l{k), l{p))d{p + k-p'- k'). (10) 

All the non-mixed species evolve according to the standard Boltzmann equation 



ni{p) = / dp'dkdk' 



J2F,Apk\p'k')[njik')nj{p') - (11) 

3 

+G,{kp\k'p')'PT{p)-P*T{k)], (12) 

where again goes over all species. As before there is a new term which arises from annihilation of mixed states. 
The terms in Eq. (^ containing c, C, D and d are damping terms from elastic and inelastic scatterings which break 
coherence of the oscillations. Details of how to calculate these terms can be found in Ref. [ p2[ . 

Finally, the V x P term is the usual oscillation term which is responsible for the flavour oscillations. The potential 
vector V can be written as 

V - 2£;„^x + {E^o. - Epp)i. (13) 

Here, 

Eaf3 = sin 260 - Vc^pip) (14) 
2p 

Ec^a - Epp = -— COS20O + Vo.{p) - Vp{p), (15) 
ip 

where 6m? — rn^ — m\ and 6q is the vacuum mixing angle. The matter potentials, V , arise from the neutrino 
interactions with the medium p3| , p^ . 

The above set of equations is quite complicated to solve, both because of its non-linearity and because it is momen- 
tum dependent. However, it is simplified enormously by using the so-called quantum rate equations instead of the full 
quantum kinetic equations presented above. In the quantum rate equations, the original quantum kinetic equations 
are integrated over momentum space so that the momentum dependence disappears. This integration can only be 
accomplished analytically if one assumes kinetic equilibrium for the neutrinos, and therefore involves an assumption. 
We assume that the one-particle distribution functions for neutrinos (the diagonal parts of the density matrix) are of 
the form e~P/'^ . 

The usual number densities of the mixed neutrinos are then related to the integrated density matrix by = 
i[Po + Pz\ and np — ^[Po — Pz\- However, this number density is not a dimensionless quantity. In order to make P 
and Po dimensionless we instead work with the dimensionless quantities P* = P/^^i^o and Pq = Po/n,yg, where n^^ 
is the number density of a decoupled neutrino species. As will be explained at the end of the section this also has the 
advantage of making the expansion of the universe simpler to treat. For simplicity we will in the remainder of the 
paper refrain from denoting P* and Pq,* with an *. 

The quantum rate equations then take the form p2] 
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where 



P = V X P - DFt - CF*^ + - Rp]z 

Pq = Ra + Rf3 



TT-j = ^ [hjUjUj - TliTli] + GiPt ■ Pt- 



(16) 
(17) 



(18) 
(19) 



hj — 1 for neutrinos and hj = j for electrons. The new parameters C, D,V, Fijand Gi are then defined as 



D 
V 

F ■ 



C{p)n{p)~ j dp'c{p,p')n{p') 
dp D{p)n{p) — J dp'd{p,p')n{p') 
dpV {p)n{p) 

dpdp dkdk' Fij {pk\p' k')n(p)n(k) 
dpdp'dkdk'G,{p'k'\pk)n{p)n{k) 



(20) 
(21) 
(22) 
(23) 
(24) 



Some details about how to perform the phase space integrals are given in the Appendix. 

In the present calculation we are interested only in a differential effect, i.e. the difference between the actual 
neutrino density and the density if neutrinos were completely decoupled. In that case, the quantum statistics of the 
involved particles is not very important |^]. In the following we therefore approximate all quantum statistics with 
Maxwell-Boltzmann statistics. This means that Bosons and Fermions have exactly the same behaviour. 

In previous calculations of neutrino oscillations in the early universe it has been assumed that the active neutrino 
species have the same temperature as the electromagnetic plasma. However, that is not the case during electron- 
positron annihilation. We therefore have to operate with different temperatures for all the different species. Specifically 
we always assume a Maxwell-Boltzmann distribution with a temperature given by Toffji = TQijii/riug)^/'^ , where n^^ 
and To are the number density and temperature of a completely decoupled non-mixed neutrino species, is the 
actual density of the given species. Electrons, positrons and photons are assumed to be in full thermal equilibrium 
with the temperature = T^. 

All of the above equations have been derived assuming a non-expanding universe. If the universe expands then for 
the usual one particle Boltzmann equation, the Liouville term df /dt is changed from df /dt to df /dt~ Hpdf /dp pl| . 
Likewise, the left-hand side of the integrated Boltzmann equation is changed from ?i to n -I- 3Hn [|l|. However, the 
equations can be made to look exactly like they do for the non-expanding case if they are recast in comoving quantities. 
The momentum of a particle redshifts with expansion as p cx R~^. For the momentum dependent Boltzmann equation 
one then defines the comoving momentum as = pR, a quantity that does not redshift. The Liouville operator then 
becomes df /dt = df /dt — Hpdf /dp = df{p<,)/dt, i.e. it looks exactly like the non-expanding case. 

For the integrated Boltzmann equation this also holds. The usual procedure is to write the number densities in 
units of entropy density as = n/s. If entropy is conserved n* is not affected by cosmic expansion. However, in 
the present case entropy is not conserved because full thermodynamic equilibrium is not maintained. Instead one 
can rescale the number density in units of a completely decoupled neutrino species, n^g cx R~^ . In this case the left 
hand side of the Boltzmann equation reads n -I- ZHn = . If one recasts the quantum rate equations in units of n^^ 
the cosmological expansion does not appear anywhere, and one can readily use Eqs. (p^p4[) even in an expanding 
universe. 



A. Matter potentials 



The diagonal part of the matter potential comes both from interactions with other species, and from self-interactions. 
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The electron lepton number is known to be small (Lg = Lp ^ 10"^'') because of charge neutrality. However 
the neutrino lepton numbers are not well constrained at present. At present the strongest constraints come from 
considering Big Bang nucleosynthesis (BBN) and Cosmic Microwave Background Radiation (CMBR) arguments. 
From BBN considerations one finds an upper bound on relativistic energy density (usually expressed as the effective 



number of neutrino species Ni, = p/ Pu,,) |25|2^, and therefore also on neutrino lepton numbers |27|. However, 
because electron neutrinos enter directly into the weak interactions that regulate the neutron to proton ratio an 
electron neutrino chemical potential cannot be directly translated into an effective number of neutrino species [ p7[ . 
CMBR on the other hand is not sensitive to neutrino flavour but only to energy density. Therefore the CMBR bound 
on the effective number of neutrino species can be directly translated into a bound on neutrino lepton numbers [p8|-^ . 
Hansen et al. recently combined BBN and CMBR to derive the tightest present constraint of —0.01 < Li,^ < 0.22 
and \L^^ ^ \ < 2.6 (see also p^-p7|]. This constraint is of course many orders of magnitude larger than the known value 
of the electron lepton number. 

For the sake of simplicity we assume that all lepton numbers are of the same order as Lg. In that case they can 
be ignored in the calculation. This simplifies the numerical computations tremendously because the neutrino and 
anti- neutrino sectors decouple (i.e. the equations describing them are identical). However, we stress that this is not 
necessarily a good approximation, depending on the actual values of neutrino chemical potentials. 

If one neglects all lepton numbers, the result for electron neutrinos is 

' ^ (25) 



y „ 



T,T^ + \{l~x)Tl 



where x = sit? 9w — 0.226. For the muon neutrino one finds a similar expression 

96V2Gi. 



w 



(26) 



This is very close to what is found using Fermi-Dirac statistics for neutrinos. For FD statistics the front factor should 
be 45^(31 — 9.68 instead of 96/71^ ~ 9.72 for MB. More details about the calculation of the matter potentials 
can be found in the Appendix. 

In addition to the diagonal part of the matter potential there is an off-diagonal part due to neutrino-neutrino and 
neutrino- antineutrino forward scattering This term is proportional to /P(p)dp = P, so that in the quantum 

rate approximation the term is identically zero. The off-diagonal term has the effect of synchronising the oscillations 
of different neutrino modes. However, in the quantum rate approximation the problem is reduced to following a single 
"effective" mode and therefore the oscillations of different modes is in some sense already synchronized. In fact there 
is an additional off-diagonal term which is proportional to pap , but as will be explained in the Appendix this term is 
always very small. 



B. Annihilation and damping terms 



The annihilation term for the process ii jj for particle i can be written as 



Rij — Fij [hj 



rijUj — n^.rijj. 



j 
1^0 



i 

7^8 

1^0 



FnB 



j 

7^8 

I'D 



rp>i 

i 

7^8 

I'D 



(27) 



where B depends on the specific process. Table II lists the value of B for different processes. Again, the front factor 
value 4 is very close to what is found using FD statistics (~ 3.97) |Q. 

The same front factor can be used to calculate the damping coefficients C and D. The calculation of these terms 
is discussed in Ref. [p2l for the case of FD statistics. Here we just state the result for MB statistics 



D 



1 

2^° 



4rp4 



7^8 



7^8 



rp4rp4 

7^8 



2(^ 



rp4rp4 



+ (8a;^ + 4a; + l)- 



7^8 
7^8 

I'D 



7^8 

rp4 rp4 

7^ 

1^0 



(28) 
(29) 
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TABLE II. Values of B for different annihilation processes. 



Process B 

v^Ve ^ e+e^ 8x^ + 4a; + 1 

l^lJ.^'^^ e+e~ — Ax + 1 

VrVr <^ e^e~ — 4a; + 1 

Z/eZ^e ^ !/^i7^ I 

VeVe ^ VtVt 1 



ji4 ji4 

C = 2(16a;2 + 4)Fo^^ (30) 

I'D 

Ge = (32x2-4)Fo^|^ (31) 
G.. - 4Fo^;^. (32) 

fa 

These expressions are very similar to those derived by McKellar and Thomson |^ who used FD statistics but assumed 
identical temperatures for all species. 



C. Time-temperature relationship 

Fundamentally, two equations are needed to fully describe the cosmological expansion with time ||^ . The system of 
photons and electrons/positrons can to an extremely good approximation be assumed to be in full thermal equilibrium 
via electromagnetic interactions, so that they can be described by a common temperature, T^. Therefore the two 
variables describing the cosmological expansion with time can be taken to be the scale factor, i?, and the photon 
temperature, T^. 

As the two independent equations we choose the Friedman equation |^ 

= ^ (33) 

and the equation of energy conservation |2j 

d{pR^) + Pd{R^) = (34) 

as our fundamental equations. These two equations can then be rewritten as equations for and R. Details of how 
to solve these equations in the case of Maxwell-Boltzmann statistics can be found in Ref . M . 



III. NUMERICAL DETAILS 

We have solved the quantum rate equations, Eqs. (p^p^, together with the expansion equations Eqs. ( p3[]3^ ), for 
the case of Ve — t^fi oscillations. Vt is in the present calculation assumed to be non-mixed. As initial conditions we 
choose — — 15 MeV, a temperature well above the electron-positron annihilation temperature Tann ~ 0.3 MeV. 
Furthermore we set = Py — Pz — and Pq = 2 (and identically for anti neutrinos). However, the outcome is not 
sensitive to initial conditions in the mixed neutrino sector because at high temperatures P is quickly driven to zero 
because of fast interactions, no matter what its initial value is. The system of equation is then straightforward to 
solve. 

Fig. 1 shows the evolution of the quantity 

Sn ^ - n^g ^^^^ 
n n„g 
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for Ve,v^ and Vr, for different values of the mixing angle, all calculated for the specific value of the mass difference 
= 3 X 10~^eV^. Again, n^fi is the number density of a standard, decoupled neutrino species. This figure also 
shows the limiting case of non-mixed neutrinos (sin 26*0 — 0). 

It should be noted that the result of the non-mixed case is quite close to the result found by elaborate momentum- 
dependent calculations. The calculation by Dodelson and Turner used the same approximation as we have in 
the present work (i.e. zero electron mass and Boltzmann statistics for all particles), except that they solved the full 
momentum dependent Boltzmann equation for the non-mixed case. Their result was approximately 6pu^/ Pvo — 0.012 
and 5pu^/ Pua — 0.005, whereas our result for the non-mixed case is ^Pv^l Pvo — 0.0124 and Sp^^/p^g ~ 0.0057. Notice 
that 6p/p = {pi, — Pva)l Pva is not equivalent to 5n/n^ the quantity shown in Fig. 1. However, since we assume thermal 
MaxwcU-Boltzmann distributions for all particles the two quantities are simply related by 5p/ p ~ ^5n/n. 

This shows that although our calculation is quite crude in the sense that it does not fully account for momentum 
dependence it yields results which are fairly close to momentum dependent calculations, at least for the case of 
non-mixed neutrinos. We do expect the same to be true for the mixed case, although this remains to be verified. 

Calculations which have used exact quantum statistics and electron mass give slightly smaller neutrino heating. 
Hannestad and Madsen found 0.0083 for and 0.0041 for [||, whereas Dolgov, Hansen and Semikoz found 0.009 
and 0.004 ^j^. The most recent treatment by Gnedin and Gnedin found 0.0097 and 0.0062 |||. 

In Fig. 2 we show the evolution of Pj — n^^ — n^^ for different values of the mixing angle. It is clear that for 
this specific choice of Sm^ a very large mixing angle is needed to achieve a noticeable effect. 6m^ = 1 x 10^^ cV^ 
is of the same order of magnitude as the best fitting solutions for the LMA {Sm^ = 4.5 x 10~^ eV^) and SMA 
{Sm^ = 4.7 X 10~^ eV^) solutions to the solar neutrino problem. Therefore it is clear that for the SMA solution where 
the best fit at present is sin 29q ~ 0.04 there will be no noticeable effect of neutrino oscillations on neutrino heating. 
On the other hand for the LMA solution sin 26*0 — 0.91 which is large to enough to give a significant change. 

The evolution of Pz with temperature can be understood as follows. At high temperatures oscillations are suppressed 
by the matter potential, and therefore Pz evolves independently of the mixing angle. However at a certain temperature 
the matter mixing angle goes through a maximum and oscillations become an important equilibration factor. The 
rate of equilibration between the two species is to a rough approximation given by Feq ~ D cos^ 29 sin^ 29 |4^ . This 
should be compared with the rate at which the abundances are driven apart by interactions with the electron-positron 
plasma, Fdrivc ^ (-Re ~ Rp,)- The matter mixing angle is given by the expression 

sin^ 261 = sin^ 26lo/(l - 2/cos20o + f), (36) 

where / = 6T(Ve — V^)/(5m^. This mixing angle goes through a maximum when / = cos2(?o which also corresponds 
to a maximum in the equilibration rate. This maximum occurs exactly at the temperature where the dip in Pz is 
seen. Below this temperature the mixing angle approaches the vacuum value and Foq /Fdrivc approaches a constant 
value, which to a reasonable approximation is 

This asymptotic value is always smaller than one so that for small temperatures Pz follows the same trend inde- 
pendently of the vacuum mixing angle because the driving term is dominant. However, as the vacuum mixing angle 
increases the equilibration around the maximum of the mixing angle becomes more and more important. Therefore 
the final value of Pz decreases strongly with increasing vacuum mixing angle. 

Oscillations in general become important once the matter potential no longer dominates the vacuum oscillation 
term. This happens when 

Sm'^ cos 2^0 ^ 



^^^^ < U.OxlO-^eV^j ^^^^ 

For the large value of Sttt? used in Fig. 2 this temperature is above the decoupling temperature for neutrinos. There- 
fore once the oscillations become important they are quickly damped and no oscillation pattern is seen at lower 
temperatures. 

This is not the case for lower values of Sm? where oscillations only become important well after neutrino decoupling. 
In Fig. 3 we show the evolution of Pz for various values of 5m^. For Sm^ = 10~^eV^ the oscillations are completely 
damped away because neutrinos have not decoupled before oscillations become important. For Sm^ — 10^^ eV^ 
oscillations are apparent, but still of very low amplitude. Notice that at low temperatures the mean of the oscillating 
curve still rises slowly due to neutrino heating by the electron-positron annihilations. Finally, for Sm^ = 10^^° eV^ 
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neutrinos have decoupled completely before oscillations become important. This means that there is effectively no 
damping of neutrino oscillations once they commence. 

Therefore, for low values of i5m^, Pz does not approach a definite value for low temperatures, even though the total 
neutrino number density, Pqi does. 

Of course the oscillation of Pz at low temperature is an artifact of the quantum rate approximation. In the 
momentum dependent treatment different modes oscillate with different frequencies and oscillations decohere because 
of this effect. However, it should be noted that if neutrino lepton numbers are significant, this picture can change 
completely In that case the off-diagonal elements in the neutrino matter potential are large, and neutrino 

oscillations become coherent. As long as all lepton numbers can be neglected, as was assumed in the present treatment, 
it will make better sense to use the average of Pz instead of the actual value because the decoherence effect was not 
accounted for. 

The increase in oscillation amplitude seen for decreasing 5m? is also seen for active-sterile oscillations, for the same 
reason [Q. 

One non-standard feature in the active-active oscillations is the appearance of the interference terms in the Boltz- 
mann equation. These terms could potentially be important and lead to a different result for neutrino heating. How- 
ever, it turns out that they are always completely negligible compared with the usual collision terms. Fig. 4 shows 
the quantity Q = \^iGiPT ■ Pt/ [^i'^j^j ~ "I'i^pJ for * = for the specific case of Sw? = 3 x 10~^eV^, 
sin 26*0 = 0.5. This clearly shows that the non-standard terms from annihilation of mixed states is tiny compared with 
the standard Boltzmann terms and that they can be safely ignored in numerical treatments. 



IV. DISCUSSION 



We have solved the Boltzmann equations governing the evolution of neutrinos around the time of their decoupling 
from equilibrium, including effects due to mixing of and j/^. We always assumed that the tau neutrino is unmixed. 

As could be expected oscillations have the effect of bringing the and abundances closer together. Strong 
oscillations also have an effect on the total neutrino energy density after electron-positron annihilation. The reason 
is that without oscillations, most of the heating is to the sector because electron neutrinos have charged current 
interactions. Oscillations drain away electron neutrinos into muon neutrinos, and therefore the back-reaction Vei^e 

e~ decreases in efficiency. The end result is a slightly larger neutrino energy density. As is customary we parametrise 
the neutrino energy density in units of the energy density of a decoupled massless neutrino, iV^ — Pu/ Pvq (so that 
in the absence of neutrino heating N^, — 3). However, it is also necessary to account for the slightly lower photon 
temperature if neutrino heating is accounted for, because all quantities should be measured relative to the actual 
photon temperature. We therefore use the definition 

N, = ^^ (39) 

Pvq P-f 

In Fig. 5 we plot the at low temperature. From this figure it can be seen that neutrino oscillations can change the 
effective number of neutrino species by about 1.2 x 10^"^. The effective number of neutrino species without oscillations 
is Ni, = 3.0467, and in the limit of large mass difference and mixing angle it approaches iV^ = 3.0479. This value is 
somewhat higher than what is found in the more thorough calculations using FD statistics and the full momentum 
dependent Boltzmann equation {N^ = 3.028 §, = 3.034 = 3.032 simply because of the larger 

neutrino heating when MB statistics is used (our value fits very well with that of Rcf. 0, who found — 3.046 
using MB statistics). 

For BBN calculations, there is an additional effect which must be considered. The electron neutrinos have a different 
effect on BBN than muon or tau neutrinos because they enter directly in the /3-reaction which regulate the neutron 
to proton ratio. An increased number of electron neutrinos and anti-neutrinos have the effect of increasing the n ~ p 
conversion efficiency. This in turn leads to a lower neutron to proton ratio at helium formation and in turn a lower 
helium abundance. This effect works in the opposite direction of a simple increase in N^. Furthermore, because of 
energy conservation the photon temperature is slightly lower if neutrino heating is included because photon heating 
by e+e~ annihilation is slightly smaller. This also has the effect of lowering the n — p conversion rate because there 
are slightly fewer electrons and positrons. 

To get a feeling for how neutrino oscillations change the helium production we have modified the Kawano nu- 
cleosynthesis code to take into account neutrino heating. We have then performed the numerical calculation 
for oscillation parameters corresponding to the best fit for the LMA solar neutrino solution (sin20o = 0.908 and 
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TABLE III. Change in helium abundance due to neutrino heating, for the case of 7710 = 5 x 10 



Solution 




sin 29 a 




5py^/pya 


^PvtIPvo 




LMA 


4.5 X 10"" 


0.908 


0.0097 


0.0087 


0.0057 


2.0 X 10"* 


No osc. 






0.0124 


0.0057 


0.0057 


1.2 X 10""* 


Maximal mixing 






0.0093 


0.0093 


0.0058 


2.2 X 10~* 



6m? /eV^ = 4.5 x 10 ^). The result is shown in Table III, together with our result for the case of no oscillations and 
the limiting case of large 5m? and sin2^?o- 

This shows that more helium is produced for the oscillating case. The main reason for this is the lower density 
of electron neutrinos in the case of mixing, but the increase in the effective number of neutrino species also leads to 
an increase in helium production. Interestingly, the change in helium abundance due to oscillations is of the same 
order of magnitude as the total effect. The reason is that in the non-oscillating case the effect on helium is very small 
because of cancellations. The increase in electron neutrino temperature is roughly compensated by the increase in 
the effective number of neutrinos, as well as the decrease in photon temperature. Oscillations destroy this accidental 
cancellation and therefore have a large effect. Our finding for the non-oscillating case fits well with other calculations. 
Using the same approximations as in the present paper Fields, Dodelson and Turner [^j found AYp = 1.1 x 10^^. 
More sophisticated methods find similar values in the range AYp = 1.1 — 1.5 x 10^^ 

In terms of energy density the changes to neutrino heating by neutrino oscillations are quite small. If indeed the 
large mixing angle solution turns out to be correct then the effective number of neutrino species is changed by roughly 
8 X 10"^ compared to the non-oscillating case. From CMBR the present bound on the effective number of neutrino 
species is Ny < 13 [^,^, i.e. more than two orders of magnitude larger than the effect induced by neutrino heating, 
and about 10"' times bigger than the change induced by oscillations. With precision data from upcoming satellite 
experiments such as the Planck Surveyor it could be possible to measure AA''^ ~ 0.04 which is comparable to the 
effect from neutrino heating. Even so the small difference induced by oscillations will likely remain undetectable. 

For BBN the change due to neutrino heating is of the order lO"**. At present the observational uncertainty on the 
primordial helium abundance is about (j{Yp) 0.005 which is about 50 times larger than the change. Here, 
however, the change due to neutrino heating is more significant, comparable in magnitude to the total neutrino heating 
effect. It is perhaps conceivable that the difference in neutrino heating from including oscillations could be detected. 

Finally, we stress that the present treatment is by no means definitive. A proper treatment of momentum dependence 
is missing, as is the inclusion of three-neutrino oscillations. Inclusion of momentum dependence is not likely to have 
a big effect for the relatively large mass difference and mixing angle characterising the LMA solution, but may 
be important for the vacuum solution. If I'^and v^- are maximally mixed, as is indicated by atmospheric neutrino 
measurements, then they will behave effectively as one species during neutrino decoupling. This should have the 
effect of lowering the electron neutrino temperature more than for two-neutrino oscillations, while increasing the 
overall effective number of neutrino species slightly. Altogether this amounts to the same effect as for two-neutrino 
oscillations, but slightly larger. 

It is simple to calculate the extreme upper limit on the oscillation effect by considering a complete coupling of v^, and 
Vr- In Table IV we show results of the same calculation as in Table III, but now assuming an infinitely tight coupling 
between muon and tau neutrinos. This is likely to be a very good approximation to the true state of affairs, because 
the preferred values of the mixing parameters for z^^ — ly^- mixing is at present |^ {Sm'^ ~ 2 x lO^'^eV^, sin20o — !)■ 
This means that oscillations become important already when T ~ 5 MeV, long before neutrinos decouple and also 
long before neutrino heating commences. Therefore in this case, z^^ and i>t should be treated as having effectively the 
same temperature. Note that this is only a good approximation when ^13 is small so that there is little direct mixing 
of Ve and Vr- Observations indeed indicate that this is the case. 

Maximal coupling between the muon and tau neutrinos therefore lead to a slightly larger increase in energy density 
due to neutrino heating. In the case of maximal ly^ — coupling it increase Ni, from 3.0479 to 3.0484, and for the 
LMA solution from 3.0476 to 3.0478. Helium production is also slightly increased, by about 0.2 x 10^*. Although we 
have not performed a full three-neutrino oscillation calculation, this estimate should be fairly close to the true value 
because and Vr are most likely maximally mixed with a large mass difference pO| . 
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TABLE IV. Change in helium abundance due to neutrino heating, for the case of 7710 = 5 x 10 . Results are from assuming 
infinitely tight coupling between i/p and Vr. 



Solution 




sin 29q 








AYp 


LMA 


4.5 X 10"'' 


0.908 


0.0089 


0.0077 


0.0077 


2.3 X 10"* 


No osc. 






0.0124 


0.0057 


0.0057 


1.2 X lO"* 


Maximal mixing 






0.0082 


0.0082 


0.0082 


2.5 X 10"* 
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V. APPENDIX 



A. Matter potentials 



The diagonal matter potentials arise from neutrino loop interactions with the background medium. The magnitude 
of the potentials were first calculated by Notzold and Raffelt Q (see also Q). In the absence of lepton numbers, 
there is no contribution due to interactions with neutrinos of different flavour. Neutrinos of the same flavour yield 
the contribution 



3m| 



(40) 



For electron neutrinos there is an additional contribution from interactions with the background electrons and 
positrons. If the electron mass is neglected the contribution is 



(41) 



Neglecting the electron mass in the matter potentials only leads to very small errors and is consistent with neglecting 
it in the interaction matrix elements. In the quantum rate approximation one should make the replacement p — *■ 
(E^,) = 3T^,. Then using = 2T^/n'^ and = T^/tt^, one finds 



96\/2G'f 



T.T^ + \{l-x)Tl 



(42) 



for the electron neutrinos. For the muon neutrino there is no contribution from electrons, and one finds the result 



(43) 



As was mentioned in Section II. a there is an additional off-diagonal term which is of the form 



V = 



3m| 



(44) 



However, because the off-diagonal elements of pap are always very small until long after neutrino freeze-out this 
contribution is negligible. 



10 



B. Phase-space integrals 



In this section we discuss how to perform the phase-space integrals needed to calculate the terms in Eqs. (P4). A 
full derivation of all the terms would be too lengthy, but as a representative example we show the calculation of the 
contribution to Re by the process Vei^e e~^e~ . 

1 



Re{VeVe 



<-> e e = 



dpdp'dkdk' {2n fV^{pk\p'k')S{p + k-p' -k') 



[fe{p')fe{k')-UAp)UAk)\ 



(45) 



Here, dp = 2po{2-nY ' ^'^^^'^ ~ ^ p /^-s, and f^Sv) — ^ ^Z^''^. All quantities are normalized to the density of a single 
decoupled neutrino species, which explains the in front. 
The squared matrix element is 



V^{pk\p'k') = ?,2Gl[{2x + lf{p ■ k'){k ■ p') + {2xf{p ■ p'){k ■ k')\. 



(46) 



Because of the Boltzmann statistics, fe{p')fe{k') — fe{p)fe{k)- Using this, the dp',dk' integrals can be performed 
using Lenard's formula 



(47) 



The result is 



1 16 ttG? 



Bx^ + 4x + 1) 



3 (27r)8 
J 2po2ko^ 



(48) 



We then use that {p ■ fc)^ = p^fc^(l — cos0)^, where 9 is the angle between the direction of p and k. After performing 
the integrals over d^p and d^k the result is then 



ReiVe^e ^ 6+6 ) 



1 4(^2 



(49) 



Using n^a = /tt^ the contribution to the Boltzmann collision integral is 



4(^2 



3 T^j8x2+4a:+l) 



7^8 7^8 
^ e _ t/e 
7^8 7^8 



(50) 



This result is almost identical to what is found using FD statistics (3.97 in the front factor instead of 4). 
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FIG. 3. The temperature evolution of Pz = n„ 
is for 5m^ = 1.0 x IQ-'^eV^, the middle for Sm'^ -- 



— n,j^ for sin26o = 0.9 for various different values of 5m^. The upper panel 
1.0 X 10"**eV^ and the lower for Sm'^ = 1.0 x lO'^'^eV^ 
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FIG. 4. The evolution of the parameter Q = | G^Pt • Pt/ X/j ~ "■"^i'^i'i] fc>r i = i/e, for the specific case of 

W = 3 X 10"'^ eV^ sin 26*0 = 0.5. 
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FIG. 5. The total change in effective number of neutrino species i^Ni, — N^, — 3) as a function of Srn^ and sin 26^0. Shown 
are also the best fit values for the possible solar neutrino solutions ll£f|. 
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